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Abstract — We predict the generation of noncritically squeezed 
light through the spontaneous rotational symmetry brealting 
occurring in a Kerr cavity. The model considers a x'^' cavity 
that is pumped by two Gaussian beams of frequencies loi and 
L02. The cavity configuration is such that two signal modes of 
equal frequency lus = (t^i +(J2)/2 are generated, these signal 
fields being first order Laguerre-Gauss modes. In this system 
a spontaneous breaking of the rotational symmetry occurs as 
the signal field corresponds to a Hermite-Gauss TEM mode. 
This symmetry breaking leads to the perfect and non-critical 
(i.e., non dependent on the parameter values) squeezing of the 
angular momentum of the output TEM mode, which is another 
TEM mode spatially orthogonal to that in which bright emission 
occurs. 

Index Terms — quantum fluctuations, four-wave mixing, non- 
linear optics, squeezed light. 



I. Introduction 

Squeezed light is a kind of radiation exhibiting reduced 
fluctuations with respect to vacuum in some special observ- 
able. This occurs at the obvious expense of an increase in the 
fluctuations of its canonical pair, as followed by the Heisen- 
berg uncertainty relation satisfied by the couple. In a single 
mode field these canonically related observables correspond 
to orthogonal field quadratures, which are equivalent to the 
position and momentum of a harmonic oscillator Squeezing 
is a macroscopic manifestation of quantum phenomena that is 
attracting continuous attention since the late seventies of the 
past century [JJ, [2 1, [3}. Nowadays a renewed interest has 
arised because of the importance of squeezing in generating 
continuous variable entanglement, which is a central issue for 
continuous variable quantum information purposes [4|. 

Squeezed light is generated by means of nonlinear optical 
processes, such as parametric down-conversion or four-wave 
mixing. The squeezing level attainable in such nonlinear 
optical processes depends on the interaction time that is limited 
by the nonlinear medium length. Thus in order to increase 
the squeezing level, these processes are usually confined to 
occur within an optical cavity. In this way the squeezing level 
can reach the largest possible levels at the system bifurcation 
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points such as, e.g., at the emission threshold. Squeezing levels 
as large as 90% (lOdB reduction respect to vacuum fluctua- 
tions) have been recently reported ||5], IS) in such conditions in 
degenerate optical parametric oscillators (DOPOs). However 
perfect squeezing (i.e., the complete suppression of quantum 
fluctuations in a field observable) cannot be achieved in these 
conditions because complete suppression of fluctuations in a 
mode quadrature implies the existence of infinite fluctuations 
in the other quadrature, what would require infinite energy in 
the process. 

Nevertheless perfect squeezing could be actually produced 
and we have recently proposed a way for obtaining it Q. 
The idea can be put in short as follows. Consider a nonlinear 
optical process in which two photons with equal frequency are 
generated, each photon corresponding to ±1 orbital angular 
momentum (OAM) Laguerre-Gauss mode. This is equivalent 
to generating two photons in a TEMio Hermite-Gauss mode 
whose orientation in the transverse x-y plane is determined 
by the phase difference between the two Laguerre-Gauss 
photons, let us denote it by (p. Now assume that (p is not 
fixed as it occurs, e.g., in a down-conversion process. This 
amounts to saying that the orientation of the Hermite-Gauss 
mode is not fixed as is the angle formed by the Hermite- 
Gauss mode with respect to the x-axis. In these conditions we 
can expect the occurrence of arbitrarily large fluctuations in 
the Hermite-Gauss mode orientation, which suggests that the 
canonical pair of (j), namely the angular momentum —id/dcf), 
could be perfectly fixed. But the angular momentum of a 
TEMio mode forming an angle with respect to the x-axis is 
another TEMio mode forming an angle </) + 7r/2 with respect 
to the X-axis. Then this mode could exhibit perfect squeezing 
in one of its field quadratures. Notice that the concept of 
bifurcation is not involved in this discussion and that the 
variable exhibiting arbitrary fluctuations is an angle. Then, 
a priori, perfect squeezing is possible in such a process as 
"infinite" fluctuations are possible in the fluctuations of 0. 

In Refs. |7|, |8| we have recently theoretically demon- 
strated the above ideas in a model of DOPO tuned to the 
first transverse family at the down-converted frequency. The 
requirement that the angle can take any possible value (i.e. 
that the phase difference between the two Laguerre-Gauss 
modes be arbitrary) is nothing but the requirement that the 
system be rotationally invariant around the optical cavity axis. 
Hence the resulting squeezing can be understood too as the 
result of the spontaneous breaking of this rotational symmetry, 
as the emitted Hermite-Gauss mode is obviously no more 
rotationally invariant. In ||9l we extended this study to DOPOs 
having different transverse families resonating at the down- 
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Figure 1.- Scheme of the system. A x''''' medium is confined within 
an optical cavity and pumped by two Gaussian beams of frequencies 
uo\ and lj2- The cavity tuning is such that two signal modes with 
frequency uj^, = (wi +0^2) /2 are generated. The two signal modes 
are degenerated in frequency but differ in the spatial mode, one (the 
other) corresponding to a Laguerre-Gauss mode with orbital angular 
momentum I = +1 (I = —1). 



converted frequency, arriving to the same conclusion: Every 
time the nonhnear process generates light which breaks the 
rotational invariance of the system, the expected perfectly 
squeezed observable is found. 

In the present paper we present a model for a x*^"^^- 
nonlinear cavity in which squeezing appears as the result of 
the rotational symmetry breaking. The interest of this new 
proposal is twofold. On one hand it allows us to demonstrate 
that rotational symmetry breaking is a robust means for 
generating squeezing in the sense that is not limited to a 
x'^^-'-nonlinear cavity such as the DOPO. On the other hand 
perfect rotational invariance could be problematic to achieve 
in x'"^^ systems because phase-matching requirements could 
imply the tilting of the nonlinear crystal thus compromising 
rotational invariance, a difficulty that disappears in a x*^^^ 
process because phase-matching occurs easily in this case. 

The type of x'^'-nonlinear cavity system we are proposing 
here is a novel one that has not been studied previously, as 
far as we know. Hence we must derive the quantum model 
(Section II) as well as study its classical emission properties 
(Section III) before addressing its quantum properties (Section 
IV). We are able to demonstrate that the proposed device effec- 
tively exhibits perfect squeezing originating in the rotational 
symmetry breaking. In Section V we resume our main results. 

II. Model 

Consider an optical cavity with spherical mirrors containing 
an isotropic x'"'^^ medium. The cavity is pumped from the 
outside with two coherent fields of frequencies uji and u!2, 
these pumping beams having a Gaussian transverse profile. 
Suppose, for simplicity, that these pumping beams have the 
frequencies and shapes corresponding to two consecutive 
longitudinal modes of the optical cavity. Then, within the 
cavity the nonlinear interaction generates, through a four- 
wave mixing (FWM) process, two other fields having the same 



frequency tOg such that lui + uj2 = 2ujs. Assume now that the 
cavity geometry and tuning is such that these two signal fields 
have the shape of first order Laguerre-Gauss modes. These 
modes carry OAM and its conservation imposes that one of 
the signal fields carries positive OAM with I = +1 while the 
other carries negative OAM with Z = — 1 as the pumping fields 
have zero OAM (see Fig. 1). 

As stated, the just described FWM process requires that 
the optical cavity modes, as well as the fields' frequencies, 
be properly chosen. An immediate choice that verifies the 
previous requirements is a confocal resonator. In this type 
of cavity the resonance frequency of longitudinal mode q 
corresponding to the transverse family f = 2p + I (p is the 
radial index) is given by ifTOl 



TTC 



1 + ./ 



(1) 



where L is the effective cavity length. In this case, the pumping 
beams can correspond to two consecutive longitudinal modes 
with / = 0. The signal modes would then correspond to 
the cavity modes with indices q and / = 1, as they verify 
2ujs — ujq,o + o-!q+i.Q — 2ujq i. Certainly, in the confocal 
resonator there are other modes with frequency uj^ (an infinite 
number indeed) having larger odd angular momenta and 
belonging to other families. However we can neglect them 
by considering that these higher order Laguerre-Gauss modes 
could have larger cavity losses (what is true for low Fresnel 
number cavities) and would consequently not be amplified. 
Once we have shown that the FWM process we propose 
could be experimentally implemented, we pass to formulate 
the mathematical model of our system. 

A. The fields 

We shall assume for simplicity that the x^^^ crystal is placed 
at the cavity's waist plane and that is thin enough as to 
perform the uniform field approximation, hence neglecting any 
dependence of the fields on the axial coordinate z. Thus we 
write the total quantum field inside the cavity, at the beam 
waist, as 

E{r,t)^Epir,t) + E,{r,t), (2a) 
Ep (r, t)= iTjAj (r, t) e"*'^^* + H.c, (2b) 

4 (r, t) = iT,A, (r, t) e-*"»* + H.c, (2c) 

where H.c. stands for Hermitian conjugate; r —r (cos 0, sin (j>) 
is the position vector in the transverse plane written in polar 
coordinates; subindices p and s denote pump and signal 
modes, respectively; = hujk/ (eonL), with fc = 1, 2, s; and 
n is the refractive index (we neglect dispersion for simplicity). 
The slowly varying amplitudes are 

i,(r,t)=a,(t)Gj(r), j = l,2, (3a) 
(r, t) - a+ (t) L+ (r) + a_ (<) L_ (r) , (3b) 

with af; (t) and aj, (t) the annihilation and creation operators 
for mode k = 1,2,+,—, which verify [am (t) , aj^ (i)] = (5„j„. 
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As for the spatial dependence in Q, they are given by IfTOl 
^ 1 



G, (r) 



J = 1,2, 



(4a) 
(4b) 



for the Gaussian and first order Laguerre-Gauss modes, re- 
spectively. In all cases wj oc l/y/uJJ, j — l,2,s. Notice 
however that in the optical domain, in which oji ^ lO^^s^^, 
one can safely take oji ~ uj2 — (^s as far as L is not very 
small, which implies that the waist is very nearly the same for 
all of the involved modes. We make this approximation that, 
although not essential, simplifies some expressions below. 

For later use we need the relation between the Laguerre- 
Gauss modes and the Hermite-Gauss modes 



V2 



(5a) 



iV2 



- = V2|L±i(r)|sin(0-a), 



(5b) 



being and iJf the Hermite-Gauss modes with an orien- 
tation a and cr + 7r/2 with respect to the x-axis, respectively. 
Thus the slowly varying amplitudes at frequency ujs can also 
be written as 



is (r, t) = dc,a (t) + &s,a (t) i?;, 



with 



1 

i 
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(6) 

(7a) 
(7b) 



the annihilation operators for the Hermite-Gauss modes. Fi- 
nally, we introduce the field quadratures of these Hermite- 
Gauss modes 



(8) 



with fleer and as. a given by Eqs. (Q. 



B. The Hamiltonian 

In the interaction picture, the system's Hamiltonian can be 
written as 

H = Hq + iJcxt + Hint ■ (9) 

Hq and Hext correspond to the modes' energies and external 
injection, respectively, and are given by 

Hq = ^ hSja^^aj, (10a) 

Hc^t = it^i {a\ — CLij + ihE2 {a\ ~ , (10b) 

with 5j = {ojcj — ^j) the cavity detuning for the mode with 
frequency ujj, being ujcj the cavity resonance closest to that 
mode. In a confocal resonator this detuning is the same for 
all the modes if the relative frequency of the pump modes is 
locked to the free spectral range of the cavity, i.e., 1^2 — ^1 — 



TTc/L. Hence, in the following we take Sj = S V j. £j are the 
pumping parameters, which are related to the experimental 
parameters by 



2ujjhL 



being T {ujj ) the transmission factor at the considered fre- 
quency and Pj the power of the pumped laser In the following 
we will assume £1 ~ £2 — £ for simplicity. Hint describes the 
nonlinear interaction and can be written as the sum of three 
contributions 



-hg (^Hsp 



in I -^-'cpm 



with 



-ffspm = a\^al + Clival + - (^a+ a+ + a'Ja": 
Hcpm = 4:a[diala2 + 2a^a+aLa_ 

+ 2 (ajai + 0,202 ) [o,\-<i+ + a^ci- 



•J2f2 I ~t2-2 



^fwm — 2 (a\ald+a- + 01020^0^ 



(11) 

(12) 
(13) 

(14) 



describing self-phase modulation, cross-phase modulation, 
and four-wave mixing, respectively. Note that this Hamilto- 
nian contains all the possible combinations of four operators 
conserving both energy and OAM. The factors multiplying the 
different terms are intuitive once one takes into account the 
two following features: (i), there is a global factor 4 in iJcpm 
and iJfwm with respect to Hspm coming from all possible 
permutations of the different operators; and (ii), if any of the 
four modes is a signal mode, a factor 1/2 appears that comes 
from the transverse modes' overlapping integral. Finally, it can 
be shown that the coupling constant is given by 
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(15) 



with w the beam waist, l^. the crystal length and x the third- 
order nonlinear susceptibility of the crystal. 

C. The quantum evolution equations 

In this subsection we apply the standard procedure to 
develop the quantum theory of a nonlinear resonator within the 
generalized P representation to our Kerr cavity model [11], 
Wi\ , |13|. The starting point is the system's master equation 
for the density operator g, which reads 

_d . _ ^ r. 

dt ^ ih 

where C is the Liouvillian superoperator describing field 
losses through the output mirror, which applied to the density 
operator reads 



H,Q 



(16) 



E 



aj , ga 



ajg, a' 



(17) 



As we are assuming that the system has perfect rotational 
symmetry around the cavity axis it follows that 7+ = 7_ = 7s. 

As usual, we use now the generalized P representation in 
order to transform the operator master equation into a partial 
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differential equation for the quasiprobability distribution P. In D. The reduced model 



this representation, to every pair of boson operators (^j , aj 
it corresponds a pair of independent stochastic amplitudes 
(oij, a^) verifying {a1j) — {oij)* ■ By using standard tech- 
niques, one finds that the Fokker-Planck equation governing 
the evolution of P reads 

1 92 



d_ 

m 



P{a]t) = 



d 



dc 



2 dajdc 



Pia-t), (18) 



where we write vector a as 



a = (ai, Q!2, a+, a_, a^, a^, at) , (19) 

and give the explicit expression of the components of both the 
drift vector A (a) and the diffusion matrix V (a.) in Appendix 
A. 

Once the Fokker-Planck equation is known, it can be con- 
verted, by applying Ito rules, into an equivalent set of stochas- 
tic first-order differential equations: The quantum Langevin 
equations. They read 



— a ~ A (a.) 



Bia)-rj{t) (20) 

where the components of t] (t) are real gaussian noises veri- 
fying 



and the noise matrix S can be obtained from 



(21a) 
(21b) 

(22) 



The equivalency between the master equation and Langevin 
equations must be understood as follows 

(: / {am, al) :) - (/ , (23) 

i.e., quantum expected values of normally ordered functions 
are equal to the stochastic averages of the same functions after 
changing boson operators (a„i,aj,j) by independent complex 
stochastic variables (am,a+). 

With this procedure, we have obtained a set of stochastic 
differential equations, Eqs. (|20] | together with Eqs. ( l64l i and 
(l66T l. ruling the evolution of the fields inside the cavity. 
However, the analysis of this model turns out to be quite 
involved. 

A simpler model retaining the essential ingredients consists 
in neglecting the temporal evolution of the pumping fields, 
or further, considering the undepletion limit for the pump. As 
follows from ( I64ab and (164b b applied to (l20l i. this will be a 
reasonable approximation whenever the gain g, Eq. (fTST i, is 
small enough as for \a±\ <C |ai.2|- In such a case one can 
safely neglect the amplitude variation of the pump modes. Of 
course, in doing that some of the dynamical richness of the 
system is lost, but the approximation simplifies very much the 
analysis of the rotational symmetry breaking which is fully 
retained by the simplified model. 

Hence, in the following we will study a reduced model in 
which the pump fields are taken to be equal (and real without 
loss of generalization), i.e., we take 



ai 



a2 = p e 



(24) 



Under the assumption that the pump fields remain constant, 
thus satisfying ( l24l i. the Fokker-Planck equation of the system 
looks like Eq. ( fTSl l, but with simpler diffusion matrix and drift 
vector. 

From the general expressions given in Appendix A it is easy 
to obtain that the diffusion matrix reads now 



V 



*5 





^.2 



2p^ 



(25a) 
(25b) 



(1?^+^ is like ^ but swapping ak and a^, and changing i 
by — i), and that the components of the drift vector are 



(7s + iS) a±+ i {a^a± + 2a|a^ + 4^^) a± 

+ 2ip^a+, (26) 



being like A^^ after swapping a^. and a^, and changing 
i by — i. 

Like before, from this Fokker-Planck equation, a set of 
Langevin equations can be obtained. Nevertheless, it will be 
useful to introduce the following change of variables and 
parameters (we note that the phase factor ip appearing in the 
rescaled fields is introduced to make the stationary solutions of 
the classical equations satisfy /?_ = see Eqs. (|33] |. which 
simplifies the quantum analysis) 




a±e 



with 



sin 2^p 
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— , K = 


^, T = 


Is 


Is 




1 


9 2p 





(27) 



(28) 



in terms of which the Langevin equations for the reduced 
model read 



2i-4> 



2p+p^ + 2p) f3i 



ipe -^'f3^+Bfj^,,(j{T), (29) 

plus the corresponding equations for (3^ which are like those 
for /3± above after swapping l3k and /3^, and changing i 
by —i. The overdot indicates derivative with respect to the 
adimensional time T, and the four components of the noise 
vector ^ satisfy properties ( |2T] ) now for the adimensional 
time T. As for the noise matrix B, it can be obtained from 
V — B ■ B as usual, but now using the diffusion matrix after 
introducing the changes dZTl i. which reads as ( 125 ab but with 



2a+a_ + 2e-2»'Ap2 
^,2 



and with like but swapping ak and a^, and 

changing i by ~i. In any case, there is no need to evaluate B 
at this moment. 

Eqs. ( |29l ) constitute the model we analyze in detail below. 
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III. Classical limit 

Before addressing the analysis of quantum fluctuations we 
need to know the classical steady states of the system as well 
as their stability properties, and for doing that we must first 
write down the classical Umit of Eqs. (|29] l. This is easily done 



by neglecting the noise terms and by making /? 



the quantum Langevin Eqs. ( |29l ). We obtain the following set 
of two complex ordinary differential equations for the signal 
classical fields amplitudes (3±: 



/3± 



1 + z(A-|/3±|2-2|/3^|^ 



2p 



(30) 



Thanks to rescaling (|27] | we can appreciate that the classical 
dynamics of the system is governed by just two parameters, 
namely the normalized detuning A and pump strength p. We 
pass now to study the stationary solutions of Eqs. ( [30] l as well 
as their stability properties. 

Eqs. (|30] | have two steady states. First there is the trivial 
steady state /3± — 0. It is easy to show that its stability is 
governed by the eigenvalues 



A± = -l± Vp2-(A-2p)^ 



(31) 



what implies that the trivial solution is stable except when the 
pump amplitude p verifies p_ < p < p+ with 



P± 



\ f2A± VA2-3 



(32) 



in which case the trivial solution becomes linearly unstable 
because Re(A+) > 0. Notice that a prerequisite for the 
destabilization of the trivial solution is that A > \/3 {p is 
positive). 

At the instability points Re(A+) ~ the system passes 
from the trivial state to the nontrivial one, which reads 



with 



M2 = i(A-2p±v^^), 



(33) 



(34) 



and where 6 is half the phase-difference between the two 
Laguerre-Gauss modes amplitudes, and is not fixed by Eqs. 
(|30] |. being hence arbitrary. In Fig. 2(a) this solution is shown 
as a function of p for three values of A. 

Notice first that there are two possible values for ji. It is easy 
to demonstrate that the solution with the minus sign in front of 
the square root (dashed lines in Fig. 2(a)) is always unstable, 
while the solution with the plus sign (continuous lines in 
Fig. 2(a)) is stable within all its domain of existence. Notice 
also that while the sum of the phases of the two Laguerre- 
Gauss modes amplitudes is fixed to zero (this is thanks to 
the change of variables (|27]|). their phase difference 6 is an 
arbitrary quantity, which reflects the rotational invariance of 
the system as this phase difference determines the orientation 
in the transverse plane of the emitted Hermite-Gauss mode, 
see Eq. ( |36] | below. 




Figure 2.- (a) Intensity of the emitted classical signal beam as 
a function of the pump strength p for the indicated values of 
the detuning A. Continuous (dashed) lines indicate linearly stable 
(unstable) solutions, (b) Domain of existence and stability of the 
solutions. The non trivial solution exists in the domain limited by 
the two continuous lines. The trivial solution is stable to the left of 
the dashed line and to the right of the right continuous line. Hence, 
the system is bistable in the region between the dashed line and the 
left continuous line, marked in the figure. 



From Eq. ( |34] | it follows that the necessary condition for the 
existence of the nontrivial solution is p > 1, and its domain 
of existence is determined by the condition 



A > 2p - Vp23Y^ 



(35) 



Note that this inequality can be satisfied only if A > -\/3, in 
agreement with our analysis of the trivial solution. In Fig. 2(b) 
we represent this domain of existence as well as the domain 
of stability of the trivial solution, see Fig. 2(b) caption. 

From all the above we see that the classical nontrivial state 
slowly-varying amplitude can be written as 



Af"-'' (r) = ^ [€-'"1+ (r) + e+'^L^ (r)] 



(a - 2p + v/p^ - 1 



-,1/2 



Hi (r) . (36) 



This solution corresponds to a Hermite-Gauss mode rotated 
an angle 6 with respect to the x-axis, being the value of this 
angle arbitrary. As already noticed, the arbitrariness of 9 is 
the consequence of the rotational symmetry of the system: 
As the pumping modes are rotationally symmetric, as the 
cavity is, the signal TEMio Hermite-Gauss mode cannot have 
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any preferred orientation and thus all orientations are equally 
likely. This property is essential for the results we present 
below for quantum fluctuations. 

Finally, note that the fact that the emission takes place in a 
TEMio mode, allows us to distinguish between two different 
modes: The (r) mode, which is the classically generated 
one, and its orthogonal mode (r), which is empty of 
photons at the classical level. In the following we shall refer 
to these modes as the bright and dark modes, respectively. 

IV. Quantum analysis 

We return to the quantum description of the system. In the 
following subsections we first derive the linearized Langevin 
equations and then solve them with the method introduced 
in m, GSI and further used in |7l, fS], r9|. Next we show 
the most outstanding quantum properties of the system: After 
proving that the bright mode is rotating randomly in the 
transverse plane, we show that the dark mode has perfect noise 
reduction in one of its quadratures. 

A. Linearization of the Langevin equations 

We linearize Eqs. ( |29l ), a valid approximation in the large 
photon number limit and thus appropriate for our purposes 
as we are going to analyze quantum fluctuations around the 
above threshold solution ( [33l ). Consequently we write 



(37a) 
(37b) 



where the phase difference 6 appears explicitly in order to 
keep the b's small. Then, by assuming that bj, ^j, and 9 are 
small quantities, we easily arrive to the following linearized 
Langevin equations 



where 



b + iNon0vo ^Ch + K.B£, 



b = col (6+,6_,6+ 6i) , 
vo = col (-1,1, 1,-1), 



(38) 

(39a) 
(39b) 



A^o is a normalization factor (see Eq. (|49] l below), £ is a matrix 
with elements 



and K, is the diagonal matrix 

IC = diag (e'^ e-'"', e-*^ e'' 



(41) 



Finally, B refers to matrix B evaluated at the stationary state 
(|33l l. Its expression can be derived from V = B-B with 



V 



(-) 




(42a) 



2/i2 +pe-2'V' 



(42b) 



After some algebra, it is easy to show that B can be written 

as 



B^ 



b'-' = 



B 



{-) 



B 





(+) 



ae 



be'''' 



(43a) 



(43b) 



where elements (a, b, c, d) are independent of 6, and satisfy 
the following relations 



+ 6^ = + = iKf?, 



(44a) 
(44b) 



B. Solving the linearized Langevin equations 

In order to solve the linearized Langevin Eqs. (|38] l we 
follow a procedure analogous to that in |7|, |8|, f9\, [T4\, 
|15|. The method consists in projecting Eqs. ( |38] | into the 
eigensystem of the linear operator C. This provides a direct 
way for the evaluation of the quantum fluctuations of the 
relevant physical quantities, as we show below. 

Operator C has not an orthonormal but a biorthonormal 
eigensystem, i.e., there is an eigensystem of operators C and 
£^ verifying 



such that 



(45a) 



(46) 



The quantitative result of the analysis of C and £Ms that their 
four eigenvalues read 

Ao =0, Ai = -2, (47a) 
A2 = A2(p,A), A3 = -2-A2. (47b) 

We see that eigenvalue Aq is always null. This means that its 
corresponding eigenvector, which is said to be a Goldstone 
mode (its expression is given below), is neutrally stable, i.e., 
that its associated variable can take any possible value. Of 
course this reflects the indeterminacy of the phase difference 
between the two Laguerre-Gauss modes, 9, or what is the 
same, the indeterminacy of the orientation of the Hermite- 
Gauss output mode in the transverse plane. Hence, the null 
eigenvalue implies that the fluctuations introduced by quantum 
noise in this orientation are not damped and, hence, that 
quantum noise will induce arbitrary rotations of the Hermite- 
Gauss output mode in the transverse plane. Thus the breaking 
of the rotational symmetry of the system introduced by the 
appearance of the Hermite-Gauss mode is, in a sense, coun- 
teracted by quantum noise by making possible any possible 
orientation. Together with Ao = 0, there is the companion 
eigenvalue Ai = —2. Its corresponding eigenvector (see 
below) is consequently maximally damped irrespective of 
the system's parameter values. As we show below, it is the 
observable associated to the eigenvector corresponding to Ai 
the one that is perfectly squeezed. 

On the other hand, the two other eigenvalues, A2 and A3, are 
complex in general; nevertheless, they reach the values and 
—2 at the bifurcation points (i.e., at the points in the parameter 
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space separating the regions where the steady state solution 
(|33] l exist or not). Consequently, their associated eigenvectors 
are those exhibiting the usual squeezing occurring only at the 
bifurcations. This squeezing is not perfect (it seems perfect 
only owed to the linearized treatment) and degrades quickly 
as the system parameters are brought apart from the bifurcation 
points. We shall not analyze this squeezing here because there 
is not any relevant new feature in it with respect to what has 
been described many times in other nonlinear optical cavities. 
Hence, in the following, we concentrate on the analysis of the 
modes associated to the first two eigenvalues that are the ones 
connected with the rotational symmetry breaking. 

It is not difficult to show that the eigenvectors of £ 
associated to Aq and Ai are 



vo = col (-1,1, 1,-1), 



(48a) 



vi = — col{e''l'°,-e"^°,e-'"l'°,-e-"t'°) , (48b) 

Nq 



with 



No = -4cos(/)o, 



a normalization factor and (po a real quantity given by 



^2 +pe~*'^ 



2(^2 - (A + z)' 



(49) 



(50) 



As for the eigenvectors of £t they read 

wo = col {e"'"',-e"'"',-e- 
wi = col (1,-1, 1,-1). 



(51a) 
(51b) 



Once these eigenvectors of the linear operator are known, 
we proceed to project quantum fluctuations onto them. We 
define projections 



Cj = w* • b, j = 0, 1. 



(52) 



Note that these projections can be easily related with the 
fluctuations of the quadratures associated to modes and 
H^, see Eqs. ([8]). In particular, by using ([8j and dSTT i it is 
easy to arrive at 



Next we project the linearized Langevin Eqs. 
multiplying them by w* on the left we get 

^(r) = -^wS/C6^(r), 
ci (r) = Aici (r) + wt/c^| (T), 



(53) 
By 

(54a) 
(54b) 



where we have taken co = (this can be done because the 
arbitrary phase d can be conveniently redefined in order to 
collect the information on this mode). We notice that although 
K, and B depend on phase ICB does not as can be checked 
from Eqs. (ITTl l and (|43] |. Hence these equations are truly 
decoupled for 6 and ci. 



In the stationary limit, i.e. for large T, the solution of the 
above equations reads 



e{T)^e (0) 



i-T 



— f rfT'l (T') , (55a) 
Jo 



r 

ci (T) = w*/C6 / dT'$, {T') e-2(T'-r)^ (55^) 

Finally we evaluate, for later purposes, the correlation 
spectrum of ci. From the formal solution (I55bl i. it is straight- 
forward to prove that the correlation function of this projection 
reads 



Ci(t) = (ci (r)ci (r + r)) = 



wllCDK-wl _2|r 



(56) 



It is not difficult to obtain that vfllCVJCwl = —4k, and then 
the spectrum of correlation Ci (r) turns out to be 



Ci (iu) = / dre'^^^Ci (r) = - 



4 + ^2 



(57) 



C. Dynamics of the bright mode 's orientation 

We first analyze the dynamics of the output pattern orienta- 
tion, governed by 9. Eq. ( 15 8 at shows that the phase diffuses 
with time what means that the orientation of the classical mode 
in which emission occurs, Eq. (l36T l. exhibits a random walk. 
Then, although the mode orientation is well defined at every 
instant, it can be understood that the pattern orientation is 
undefined as after some time any value between and 27r 
could be found. This is what we understand when we say that 
the orientation of the output pattern is undetermined. 

It is important to see how much does diffuse. From Eq. 
( I55ab it is straightforward to show that the variance of 6 is 
given by 



(58a) 



59 {TY ) = dgT, 



where we have used the notation SA = A — (A), and 

" ^ 272 



= K- 



jj? sin 2(^0 + P sin [2 (0o + ^)] 



(59) 



Ap? cos^ 00 

with ?/>, ^1, and 4>q given by Eqs. ( |28] l, ( [34] i and ( |50l ). 

In Fig. 3 we represent D = de / k as a function of the pump 
strength p for several values of A. Notice that for /i ^ (i.e., 
at the supercritical bifurcation that occurs in the upper branch 
of the domain of existence of ( |33] |. see Fig. 2), D ^ 00. This is 
an intuitive result because when the output mode mean photon 
number is close to zero, the pattern orientation can be abruptly 
changed with the addition of a single couple of photons. As the 
system is brought apart from this bifurcation the mean number 
of photons rapidly increases and, consequently, it is more 
difficult for the fluctuation to change the orientation of the 
pattern, what is obviously consistent with the rapid decrease of 
D. As for the quantitative value notice that, except very close 
to the supercritical bifurcations, D is a quantity of order one 
what means that the diffusion constant dg is basically of the 
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30 



D 



20 - 



" A = 2 








' A = 3/ ' / 


- 




! / ! ^ = ^ 






H ! 1 





Figure 3.- Dependence of the normalized diffusion constant D = 
dg/hi as a function of the pump strength p for the values of the 
detuning indicated in the figure. The vertical dashed lines indicate 
the values of p where the nontrivial solution ceases to exist. 



same order of magnitude as n = gHs^ which we have assumed 
to be a very small number. Consequently, although the TEM 
output mode is randomly rotating in the transverse plane, the 
rotation is very slow except very close to the supercritical 
bifurcation points. 

D. Non-critical squeezing properties of the dark mode 

Now we focus on the main result of the present article: We 
show below that the dark mode has complete noise reduction 
on its phase quadrature irrespective of the system parameters. 
To this aim we evaluate its quadrature fluctuations as measured 
in an homodyne detection experiment. As it is well known 
(see, e.g., I16II ). quadrature fluctuations outside the cavity are 
given by the noise spectrum, which for a general problem 
quadrature (m refers to any of the transverse modes of 
our system) is given by 



(60) 



with 5",^^ {u}) the squeezing spectrum that, after taking into 
account our rescaling dZTb . can be written as 



S'^ (^) = 



+ 00 



dre 



9 J- 



{6Xl{T)6Xl{T + T)). (61) 



Defined in this way, (Cj) = means complete absence 
of fluctuations at = cj, while ((D) — \ means that 
fluctuations at = are those corresponding to the vacuum 
state. 

In the homodyning, the spatial profile of the local oscillator 
field (LOF) selects the transverse mode to be measured, 
while its phase selects a particular mode's quadrature. In 
what follows we suppose that the LOF is perfectly matched 
to the dark mode's profile i/f (r). As we saw in ( |53]) . the 
independent quadratures of this mode are xf"g and X^^^ , from 
which we can build a general quadrature as 



Xt 



1 



cos (j}Q 



Xf"g COS (fi + Ys^g sin {ip - ^iq) 



(62) 



This expression is readily obtained from the more common 
expressions X'^g — X^ g cos ip + X^Jg sin p and = 



X° g cos (?!)o 



XJg sin0o. 



Hence, by using the relation between the independent 
quadratures and the projections Cj ( l53T l. and remembering Eq. 
(ISTI i and that cq = 0, it is trivial to find the noise spectrum of 
this general quadrature, which reads 



cos^ if + sin^ (iy9 — 0o) 

cos^ 00 
sin^ (y) - ^o) 1 
cosVo l + (cj/2)^' 



(63) 



This expression shows that the quantum properties of the 
dark mode of the current system are exactly the same as 
that found in Q, El for the case of a DOPO cavity: at 

= 0, 1/°"* (^w = 0; cos^ (^/cos^^o and thus it 

has complete absence of fluctuations on its phase quadrature 
(93 = 7r/2), while another of its quadratures (ip — (j)Q in 
our case) carries only with vacuum fluctuations. Any other 
quadrature having ip between tt/2 and (po is squeezed below 
the vacuum level, though the squeezing level is smaller as p 
approaches to (pQ. These results are independent of the system 
parameters, what we expected as the noise reduction rehes on 
the rotational symmetry breaking only. 

It could seem that this result violates the uncertainty princi- 
ple, as the product of the noise spectra corresponding to two 
orthogonal quadratures is below unity. However, in |[8|, we 
have proven that this is actually not the case, as the canonical 
pair of the squeezed quadratures is not another quadrature, 
but the orientation of the dark mode 6, which is indeed 
undetermined in the long time term. 

V. Conclusions 

We have proposed a model for a Kerr cavity in which a 
spontaneous rotational symmetry breaking occurs when the 
system is beyond the emission threshold: The nonlinear cavity 
has a perfect rotational symmetry and is pumped by Gaussian 
beams, but the emitted signal field has the shape of a TEMig 
mode that breaks the rotational symmetry. We have demon- 
strated in a special simple limit (in which the pumping fields 
are taken as constants) that the rotational symmetry breaking 
implies (i) the diffusion of the output mode orientation, and (ii) 
the perfect squeezing of the phase quadrature of the TEMio 
mode that is rotated tt/2 with respect to the signal TEMio 
mode. These results are in perfect agreement with our previous 
proposal of the symmetry breaking mediated squeezing in 
a DOPO model I?], JSj. The interest of the results here 
presented are twofold. On one hand, we are proposing a system 
different to that of Q, |[8| for the possible observation of 
the phenomenon, thus showing that the results in Q, ||8] are 
quite general. On the other hand, rotational symmetry could 
be broken in a DOPO cavity if angular phase matching is 
necessary, as in this case the nonlinear crystal axis is rotated 
a certain angle with respect to the cavity axis, a problem that 
does not exist in the case of a x^'^) nonlinear medium as phase 
matching is easier to obtain. 
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As for the particular model we have proposed, in its formu- 
lation we have assumed a confocal cavity, as in this cavity type 
the required resonances are verified. Notice however that the 
ingredients that are essential for the phenomenon of rotational 
symmetry breaking mediated squeezing generation are: (i) 
rotational invariance, and (ii) that the signal field photons have 
non-null 0AM. Then, the dynamics of the pumping modes is 
irrelevant except for the quantitative details. In this sense, the 
use of a confocal cavity is not essential and any other x''^'' 
cavity in which the signal modes are the right ones could 
exhibit the described phenomenon. Then this phenomenon 
could possibly be observed in other x'"^^ cavities such as, e.g., 
fiber ring resonators. 

Another comment we would like to add concerns the non- 
pump depletion approximation. How would the dynamics of 
the pumping modes affect the results we have derived? The 
inclusion of the pumping fields equations in the study would 
obviously introduce more eigenvalues (the ones corresponding 
to the stability of these modes) and would modify some of 
the eigenvalues governing the dynamics of the signal modes, 
but would not modify the existence of a Goldstone mode 
once the signal field is switched on (as it appears due to 
the symmetry breaking). Hence as far as the dynamics of the 
pumping modes does not destroy completely the stability of 
the cw signal field emission, wherever the signal modes are 
stable the phenomenon will be present. It could well happen 
that the general model exhibits Hopf bifurcations that would 
reduce the domain of stable existence of the cw signal modes, 
but there will be a finite domain of stability for these modes 
and within it there will be the perfect squeezing properties we 
have described in our work. 

We would finally stress two important features. In |7| we 
demonstrated for a DOPO model that small imperfections in 
the rotational symmetry do not lead but to a small degradation 
of the squeezing level. On the other hand, in |8| we have 
numerically demonstrated for the same DOPO model that the 
rotational symmetry breaking mediated squeezing is perfect 
beyond the linear approximation. These conclusions should 
also hold for the system here presented. 

Appendix 

The components of the drift vector in Eq. ( fTSl l read 

A-ai — £p ~ (71 + i5) ai + 'iiga.2 a2ai+ (64a) 

Aa^ = £p- (72 + iS) a2 + 4:iga^aia2+ (64b) 
2ig {oi\oi^ + OL\a^a2 + 0^0.-0,2 + aXa+ot-'^ . 



(64c) 



2ig (a_a_Q;+ + a{\aia+ + Oi\a20ij^ + a_a\a2) , 



Aa_ = — (7s + a_ + iga^o?_^ 



(64d) 



2ig (a^a+a_ + a\a\a.^ + a^Q;2a_ + a^a.\a2) , 
and the rest of components A + are as A^^, after complex- 
conjugating and swapping and a^. As for the diffusion 
matrix in Eq. ( fTSl ), it reads 



V = 2ig 



P^-) 
-V^^ 



being ^ a 4x4 matrix with elements 

„2 



'-^11 



V 
V 



*2i -^33 



(-) _ 



= 2^32 ^ = a2a+, 2^24 = ^42 ' = "2a- 



V 



(-) 



(66a) 

(66b) 
(66c) 
(66d) 
(66e) 



"34 —-^43 — aia2 + a+a-, 
and 1?^+^ as P^^^ after swapping Uk and a^. 
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